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Abstract — We propose a contention-based random-access pro- 
tocol, designed for wireless networks where the number of users 
is not a priori known. The protocol operates in rounds divided 
into equal-duration slots, performing at the same time estimation 
of the number of users and resolution of their transmissions. The 
users independently access the wireless link on a slot basis with 
a predefined probability, resulting in a distribution of user trans- 
missions over slots, based on which the estimation and contention 
resolution are performed. Specifically, the contention resolution 
is performed using successive interference cancellation which, 
coupled with the use of the optimized access probabilities, enables 
throughputs that are substantially higher than the traditional 
slotted ALOHA-like protocols. The key feature of the proposed 
protocol is that the round durations are not a priori set and 
they are terminated when the estimation/contention-resolution 
performance reach the satisfactory levels. 

I. Introduction 

Since its inception, slotted ALOHA (SA) IT] (and its many 
variations) has been a popular distributed random access 
scheme, adopted in many wireless networks. The common 
feature for all SA variants is that, in order to maximize the 
throughput, the knowledge of the expected number of user 
transmissions per slot is required. In highly dynamical sys- 
tems, as is the typical case for wireless access, this parameter 
is usually a priori not known and it varies over time, hence 
its estimate has to be obtained for optimal operation. To this 
end, a variety of approaches have been proposed, like [[2l-||6|, 
to name a few. 

Building up on the results reported in Q, we propose a 
SA-based random access protocol suited for wireless networks 
where the number of users is not a priori known. The protocol 
operates in rounds that are divided into equal-duration slots. 
In every round users access the wireless link on a slot basis 
with a predefined probability, resulting in a distribution of 
user transmissions over slots that is used both for estima- 
tion and contention-resolution. Round lengths, measured in 
number of slots, are not a priori set - they are run until 
estimation (in the first round) or contention-resolution (in the 
following rounds) performance reach satisfactory levels. Using 
the combination of the estimation and contention-resolution 
approaches outlined in the text, we demonstrate that both 
precise estimates and high throughputs can be achieved. A 
suitable application scenario of the proposed scheme is for 
wireless access networks with large, but unknown number 
of devices, as typically encountered in Machine-to-Machine 
(M2M) scenarios. The complexity burden of the scheme is on 



the receiving side, i.e., a Base Station (BS) in M2M scenario, 
further making it suitable for M2M communications. 

The organization of the rest of the paper is as follows. 
Section introduces the background and related work. The 
system model and the principles of the proposed protocol are 
elaborated in Section |III] Section |IV] presents the simulation 
results, while Section IVl concludes the paper. 

II. Background and related work 

Framed SA in its basic variant is a simple and straight- 
forward scheme. Link-time is organized into equal-duration 
frames, which are divided into equal-duration slots. The total 
number of transmissions per frame divided by the number 
of slots, denoted as the load G, is assumed to be constant. 
Prior to start of a frame, every user randomly and indepen- 
dently chooses a single slot within that frame in which its 
transmission is going to take place. Slots can contain no user 
transmission (i.e., idle slots), just one transmission (singleton 
slots), or several interfering transmissions (collision slots), and 
only singleton slots are considered usable. The well known 
result for the throughput of framed SA is T = Ge^'^ , which 
attains the maximum T = 1/e w 0.37 when G — 1. To 
achieve G = 1 and thus maximize the throughput, the number 
of slots within the frame should be equal to the total number 
of transmissions, introducing the need for the estimation of 
the number of user transmissions. 

Far-reaching changes to the original framed SA were intro- 
duced in 1 8 1, where users were allowed to perform multiple 
transmissions (i.e., repetitions) of the same packet within a 
frame and each transmission carries a pointer to the slots where 
the other repetitions took place. On the receiving side, initially 
all transmission that occurred in a singleton slots are resolved. 
Subsequently, successive interference cancellation is executed 
based on the information carried in the pointers, as illustrated 
in Fig. [T] By this approach the throughput of the scheme 
increases; for the simple scenario when each user performs 
two repetitions in randomly selected slots of the frame [81, the 
throughput increases to T « 0.55 for logical loacQ G ~ 0.6. 

In l9|| it was noted that the execution of SIC within the 
framed SA framework resembles the execution of the iterative 
belief-propagation (BP) decoding on erasure channel, enabling 

'We note that in SA with repetitions one should distinguish between 
the logical and physical load, where the former includes only "unique" 
transmissions and the latter includes includes all repetitions. 
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Fig. 1. Successive interference cancellation in slotted ALOHA. The trans- 
mission of user 2 is resolved in singleton slot 4, enabling the removal of its 
replica from slot 1 and subsequent resolution of the transmission of user 3 in 
slot 1 . In the same way, resolution of the transmission of user 3 enables the 
removal of its replica from slot 3, thus resolving the transmission of user 1. 



the application of theory and tools of codes-on-graphs. Follow- 
ing this insight, the author analyzed the convergence of the 
SIC using and-or tree arguments |10| and obtained optimal 
repetition strategies in terms of maximizing throughput of the 
scheme. It was shown that the optimal repetition strategies 
follow the same guidelines used for encoding of left-irregular 
LDPC codes. In the asymptotic case when the number of 
users tends to infinity, both (logical) G and T tend to 1. 
Nevertheless, for the optimal performance, the number of slots 
in the frame is determined by the number of transmissions. 

Another important class of codes with advantageous 
erasure-correcting properties are rateless codes ifTTl . lfT2l . 
Rateless codes do not have a priori fixed code rate - the 
transmitter produces encoded symbols until a feedback signals 
that the message is decoded. Using the ideas analogous to 
rateless coding, aframeless SA scheme was introduced in fJ]; 
we briefly review it in the following subsection. 

A. Frameless ALOHA 

The length of the contention round in frameless ALOHA is 
not a priori fixed and the round lasts until the BS decides to end 
it by transmitting a notification to the contending users. The 
users are synchronized both on the round and the slot basis, 
implying that the numbering sequence of the slots within the 
round is a common knowledge across users. In each slot of 
the contention round users attempt to access with a predefined 
probability, termed slot access probability, which is the same 
for all users and for a given slot. The slot access probability 
for slot Sm is Q: 



Pn 



N 



(1) 



where is the expected number of interfering transmissions 
in slot Sm, denoted also as the expected slot degree, and N is 
the number of users. In the general case, /3m is a function of 
the slot number m. As shown in |7|, the above access method 
results in a LT-like distribution 1 1 1 1 of user transmissions over 
the slots of the round. This distribution can be optimized using 
standard tools from the codes-on-graphs theory. Coupled with 
SIC, for the realistic number of users the proposed access 
method enables throughputs higher the ones reported in |[9l- 

As shown in 111] and elaborated in |7|, in order to resolve 
all user transmissions with a probability that tends to 1 when 
the number of slots tends to N, the average slot degree should 
scale as 0(log N). This requirement cannot be met in practice 



even for moderately high values of N, as it would adversely 
affect the performance of the interference cancellation algo- 
rithm. Therefore, a contention round in frameless ALOHA 
ends when a sufficiently high fraction of user transmissions 
has been resolved, which can be achieved with a relatively 
low average slot degrees. The feedback that signals the end 
of the round is sent by the BS, however, unlike in the framed 
ALOHA case, the time instant at which the feedback arrives is 
a priori not known and it adapts to the contention process. The 
users with unresolved transmissions continue to contend in the 
following rounds, in the same way as in other SA schemes. 

B. Estimating the number of users 

The optimal operation of any variant of SA requires the 
knowledge of the number of users N . In many practical 
situations, TV is a priori unknown and/or it varies over time, 
thus its estimate has to be obtained. The estimation is regularly 
based on the observation of how many idle, singleton and/or 
collision slots have occurred within a frame. 

In an early work |j2l set in the framework of framed SA, 
the estimate of the current number of backlogged users is 
obtained using the number of collision slots in the current 
frame. The estimate determines the length of the next frame, 
with an aim of throughput maximization. An approach suitable 
for the cases in which the frame length is fixed is suggested 
in m, where the estimate of the number of contending users 
is based on the number of singleton slots within the frame. 
The estimate determines the probability that users choose to 
contend in/skip the next frame such that the expected number 
of contending users is tuned to the frame length and thus the 
expected throughput is maximized. 

The estimation of the number of users is in particularly of 
interest in RFID systems fSl, fSl, ["61. In f3l a simple estimator 
is derived that combines the number of singleton and collision 
slots. In [5 1 the estimation based either on the number of 
idle or collision slots, depending on the frame length and the 
number of RFID tags. The algorithm presented in |6 | improves 
the approach of |2|, enhancing the estimation by using the 
number of the observed singleton slots as well. 

All the mentioned estimators assume that the round length 
is set before it commences and that the users attempt (at most) 
a single transmission within it; this approach is unfitting for 
the frameless ALOHA scenario. Instead, we use an approach 
in which the round length is not a priori set and users proba- 
bilistically attempt to transmit on a slot basis with predefined 
slot access probabilities, as elaborated in Section UlI-AI 

III. System model 

We consider the following setup. The network consists 
of N users, where N is lower- and upper-bounded but a 
priori unknown, which contend to access the BS with their 
uplink transmissions. We focus on the resolution of a batch 
user arrival, however, the presented scheme can be easily be 
extended to the scenarios with continuous user arrivals; we 
outline the possibilities for this in Section [V] 
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Fig. 2. (a) Rounds of proposed protocol, (b) Variation of slot access 
probability within and over rounds. 



The access protocol operates in rounds that are divided into 
equal-duration slots, as depicted in Fig. |2a). The beginning 
and the end of each round is indicated by a downlink beacon, 
sent by the BS. The beacons are also used to synchronize 
the users both on the round and slot basis, to announce the 
function that determines the behavior of the slot access proba- 
bility in the next round, as well as to acknowledge previously 
resolved users. All rounds are both used for the estimation of 
the number of users N and resolution of their transmissions; 
however, the emphasis in the initial round is on obtaining 
a sufficiently reliable estimation of N. In the subsequent 
rounds the emphasis is on the contention resolution, while 
the estimation is being progressively improved. 

Let N denote the current estimate of N. Let Sij we denote 



j-th slot of i-th round. 



the corresponding slot degree (i.e.. 



the number of colliding user transmissions in slot Sy ) and 
Pij the slot access probability for Sij. The total number of 
resolved users after the i-th round (which becomes known 
as the protocol operation progresses) is denoted by Nj^f, the 
estimate of the number of contending users in the i-th round 
is thus Nci = N — N]i,i^i. By default, we assume N^q = 0. 

A. Estimation algorithm 

The BS observes each slot as being either an idle, a 
singleton or a collision slot. The probability mass function of 
this observation, conditioned on the number of users N, is 
given by: 



~riciPij{l - Pij)""^''^ if 



= 0, 

= ^' (2) 
> 1. 

and nci = 



where pij is the slot access probability for sjj 
n — nB.,i-i is the number of contending usersO 

We perform maximum likelihood estimation (MLE) of N 
using sequence of observations {s^ }. Due to the independence 
of slots, the MLE of N is given by: 

TV = argmaxTT/(sijln) = argmax ^"111/(5^ |n), (3) 



^For the sake of clarity, we note that variables n, na, riji i_i correspond 
to the random variables A'^, Nci and Nn^i-i, respectively. 



and it is obtained by solving the following equation for n: 
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Combining (|4|i and (|5]l yields: 
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where O^, Oi and Oc are sets of observed idle, singleton and 
collision slots, respectively. Eq. (|6]l is solved for n using the 
fact that nci = rt — ?ii?,i-0 this can be done efficiently using 
fast numerical root-finding methods like Brent's method [13l. 

As evaluation of ^ requires an iteration over all observed 
slots, and as the numerical root-finding methods also implies a 
separate iterative procedure, the overall evaluation complexity 
can be high. However, it is reasonable to assume that the BS 
has enough processing power to address this challenge. 

B. Protocol operation in the initial round 

The slot access probability in the first round is given as: 



Po 



(7) 



where pq < 1 and a > 1 are suitably chosen constants. The 
motivation behind Q is to progressively decrease slot access 
probability until very low values are reached, see Fig. |2|b), 
and in such way cover the interval in which the slot states will 
transient from collision to idle. The information on the number 
of contending users N is "contained" in this transient interval, 
as collision/idle slots alone can not provide information on 
which the estimate will be made |5|. 

The initial slot access probability pn = po should be chosen 
such that long (and uninformative) sequences of collision slots 
are avoided in the beginning of the initial round. On the other 
hand, the choice a determines how fast the transient interval is 
traversed, reaching the point where only idle slots appears, and 
hence N can not be further refined (implying that the round 
should be terminated). The parameter a should be chosen such 
that a compromise between a fast and a reliable estimation is 
achieved - large a implies a short round length, but also a 
short length of the transient interval, which adversely affects 

^The solution domain for n is a priori set to n > 



the reliability of the estimation. We note that the actual choice 
of pq and a depends on the expected range of the number of 
users N; these parameters are set before the initial round starts, 
and, since there is no feedback until the end of the round, they 
can not be adjusted on the way. 

The initial round should be terminated when only idle slots 
start to appear, bringing no improvement to the estimation. A 
simple rule is to end the round after observing K consecutive 
idle slots, where K is suitably chosen constant. 

The estimate N is refined after the every slot of the initial 
round, i.e., N is updated using ^ after each newly observed 
slot. Also, during the course of the initial round, some of the 
users (i.e., user transmissions) will already be resolved using 
SIC, as described in Section The output of the initial round 
is Nc2 — N — Njii, where Nm is the number of resolved 
users in the initial round. 

C. Protocol operation in the subsequent rounds 

After the initial estimate on the number of contending users 
is obtained, the main purpose of the subsequent rounds is to 
resolve user transmissions. As shown in Q, a simple and 
efficient approach is to set the slot access probability to be 
constant within the round: 

mc^) 



Pi 



Pi 



N, 



(8) 



a 



for every round i > 2, where l3{Nci) is chosen such that 
throughput of the access scheme is maximized for Nci El- 

As outlined in Section III-AI a round is terminated when 
a predefined fraction ^{Nci) of the estimated number of 
contending users Nci is resolved. In |7|, the fraction ^{Nci) 
is chosen such that the expected maximum throughput is 
reached. However, here we have to include the effect of the 
estimation errors when deciding on j{Nc'i)\ particularly as 
the choice of too high ^{Nci) coupled with an overestimation 
could lead to deadlock situations. We outline the principles for 
selecting appropriate j{Nci) in more details in Section HV-BI 

An exception to the above stopping criterion has to be made 
in the final round, when all the remaining contending users 
should be eventually resolved. A round should be designated 
as the final one when the number of contending users becomes 
sufficiently low; more details on this are given in Section lTV-BI 

IV. Results 

In this section we present simulation results for the proposed 
protocol. We assume that the number of users N is within 
range [100, 10000] and we tune the parameters of the protocol 
accordingly. We have performed only MAC-layer simulations, 
abstracting the physical layer issues (propagation, modulation, 
coding, synchronization, possible capture effects...). The SIC 
algorithm was implemented using standard iterative BP de- 
coder 1 11]. However, as outlined in |9|, despite the simpli- 
fications of the approach, the obtained results closely match 
the ones obtained by the more comprehensive and complex 
simulations, at least for moderate to low packet-loss rates. All 
results are obtained using 5000 simulation repeats for each set 
of the parameters' values. 
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Fig. 3. Normalized RMSE of the estimate after the initial round for different 
values of a, when po = 0.047 and K = 6. 



A. Choosing the parameters of the initial round 

Parameter po should be chosen such that the collision slots 
appear at the beginning of the initial round with a sufficiently 
high probability Pc, with respect to the minimal number of 
contending users A'mi„. For the targeted range of N, where 
^min = 100, Po should be chosen such that: 



l-(l-Po) 



A^minPo(l -Po) 



> Pc 



(9) 



Simulations showed that setting Pc > 0.95 does not affect the 
accuracy of the estimation; for Pq = 0.95 and A'min — 100 
we get Po « 0.047. Regarding the choice of the number of 
observed consecutive idle slots K after which the initial round 
should be terminated, the simulations showed that for K > 6 
the estimate N can not be refined any more, as measured by 
its root-mean square error (RMSE). In other words, the choice 
K = 6 seems to be a near-to-optimal one. 

Further, the simulations showed that the proposed estimator 
is not biased and that the probability density function (pdf) 
of N closely matches a Gaussian pdf that is centered at N 
and has standard deviation equal to the simulated RMSlfl 
Fig. |3] shows normalized RMSE of N, as function of N and 
for different values of a. As it can be observed, normalized 
RMSE does not change significantly with N. Fig. |4] shows the 
average number of slots of the initial round Mi normalized 
with N, for the same values of N and a. For low N the 
ratio Mi/N is rather high, which adversely affects the total 
throughput (see Section HV-Cb . Comparing Figs. [3] and |4] the 
trade-off between the accuracy of the estimation and the initial 
round length becomes clear. In further text we demonstrate the 
performance of the proposed protocol for a = 1.02. 

B. Choosing the parameters of the subsequent rounds 

The parameters of the subsequent rounds (3{Nci) and 
li^Ci) are selected using simple rules outlined below. We 
set /3{Nci) — 2.9, as this value represents a simple but 

^We leave the analytical proof of this result for future work. 
'The same holds for all the related estimators - to estimate low number of 
users with high confidence requires disproportionally high number of slots. 
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Fig. 4. The length of the initial round (in slots), normalized by the number 
of users A', for different values of ol and when po = 0.047 and /f = 6. 



Fig. 5. The normalized RMSE as function of the number of the rounds, 
when a = 1.02, po = 0.047 and K = 6; the estimate is updated only after 
each round and no backtrack is used. 



effective approximation that is based on the results given in 
ITTI . Besides providing for close-to-optimal performance of the 
collision resolution through the SIC, this choice of j5{Nci) has 
advantageous properties for the proposed estimation algorithm 
as it provides a useful "mixture" of idle, singleton and collision 
slots in the subsequent rounds (see Section Ull-Ab . 

We set ^{Nci) — 0.8 based on the following reasoning. For 
a = 1.02, after the initial round RMSE « 0.09iV (see Fig. O 
and thus N < 1.25iV and Nci < l-25^ci with confidence 
greater than 0.99. This impHes that '^{Nci)Nci = O-SNci < 
Nci in more than 99% of the cases, ensuring that the round 
is not indefinitely protracted by waiting for the non-existing 
(i.e., overestimated) users to become resolved. To address the 
remaining unlikely events of overestimation, we introduce an 
auxiliary stopping criterion - the round is forced to end after 
2Nci slots have elapsed. 

We note that setting 'y{Nci) — 0.8 does not necessarily 
imply that only up to O.SNci users are resolved within i-th 
round. Namely, it is well-known that the iterative BP decoder 
(i.e., SIC in our framework) for properly designed degree 
distribution exhibits a threshold-like behavior fT?!; in the 
beginning the number of resolved users increases slowly as the 
round progresses, however, at a certain point there is a surge 
in the number of resolved users and the fraction of resolved 
users becomes rather close to 1. In other words, in a matter 
of a single slot the resolved fraction of Nci substantially 
surpasses 0.80 provided that the actual slot access probability 
P{Nci)/ Nci is sufficiently close to the desired f3{Nci)/Nci- 
This is equivalent to the condition that Nci is sufficiently close 
to Nci', the performed simulations for the chosen po and a 
showed that this condition is satisfied with high probability 
already after the initial round and it progressively improves in 
the subsequent rounds, as displayed in Section IIV-CI 

A round i is designated as the final one when after its 
completion both of the following two conditions are met; (1) 
Pi > 0.5 and (2) \N — Nri\ < 1. The first condition ensures 
that there are sufficiently many rounds before the final one and 

^Which also denotes that this slot is going to be the last slot of the round. 



thus the estimate becomes very reliable (see Section llV-CI )n 

C. Performance of the proposed protocol 

With respect to updating of the initial estimate, we assess 
the performance of the following protocol variants: 

• iV is updated after every round i, i > 2. 

» N is updated both when number of the slots in the every 
round i reaches O.bNci and after every round i, i > 2. 

• iV is updated after every slot. 

Also, with respect to the execution of the SIC algorithm, we 
examined the following two cases: 

• The SIC is executed only on the slots that occur in the 
current round (further denoted as "no backtrack"). 

« The SIC is executed on all previously occurred slots 
(denoted as "backtrack"). 
We note that in both of the variants, we executed the SIC only 
on slots s for which |,s| < 10 (i.e., only the slots that contain 
ten or less interfering transmissions are considered exploitable 
for the SIC), as it is reasonable to assume that the interference 
cancellation can not successfully operate on high-degree slots. 
In this way the usage of high-degree slots at the beginning of 
the initial round for the SIC is avoided and the obtained results 
can be considered as more realistic. 

We are interested both in the overall estimation and 
collision-resolution performances. The estimation performance 
is expressed through root-mean-square error (RMSE) of N. 
The collision-resolution performance is measured using av- 
erage throughput of the scheme T after all users have been 
resolved; T — N/M, where M is the average of the total 
number of slots needed to resolve N users. 

Fig. |5] shows the normalized RMSE as function of the 
number of the rounds, for the case when N is updated only 
after each round and no backtrack. This could be considered as 
the worst case performance, as the estimate gets updated least 
frequently. It can be observed from that as the number of round 

^The simulations showed that the probabihty of not resolving all users when 
using the proposed criteria is minuscule. Nevertheless, an extra round/rounds 
can be appended with an aim of verifying that all users have been resolved. 
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Fig. 6. Throughput of the proposed protocol when a = 1.02, po = 0.047 
and K = 6. 

increases, RMSE drops to rather small values, implying that N 
becomes very accurate. The simulations also demonstrated that 
all the other variants of the estimation updating/SIC execution 
have only modestly better estimation performance. 

Fig. |6] shows the throughput for all the combinations of 
estimation updating and SIC execution. Obviously, using 
backtrack leads to notably higher throughputs. Further, the 
variants in which N is updated more frequently show im- 
proved performance; however, there is no difference between 
the case when the estimate is updated after each slot and the 
case when it is updated just once during the round, when 
number of elapsed slots reach half of the estimated number of 
currently contending users. As the former approach introduces 
significantly higher computation burden, it is fair to conclude 
that the latter approach is a better choice. Comparing the 
above results with the upper bound derived using the results 
presented in Q, it is clear that as N increases, the obtained 
throughputs tend to their limits. 

V. Discussion and conclusions 

A comprehensive approach to optimization of the param- 
eters of the proposed scheme, i.e., po, a, K, l3{Nc) and 
j{Nc), should take all of them into account jointly with 
an aim of the total throughput maximization. However, the 
corresponding analysis would inevitably be intractable, due to 
the multidimensionality of the problem and mutual interdepen- 
dencies of the parameters. The simplified approach used in the 
paper demonstrated the potential of the proposed scheme and 
obtained results comparable to the upper bounds given in fT\; 
a more complex optimization could be worthwhile for lower 
number of users. Nevertheless, one should be aware that the 
throughput losses for low N are mainly due to the length of 
the round where the initial estimate has to be obtained, which 
would be the case for any estimation algorithm of the type. 

When comparing the results for the variants with back- 
track/no backtrack, the former show improved performance. 
On the other hand, the memory requirements in both cases 
are comparable, as the length of the second round (when most 



of the users become resolved) constitutes the highest portion 
of the overall length. 

The integral part of the scheme are the beacons sent by the 
BS with a purpose of ending the current round/starting the 
next round with the updated slot access probability. In FDD 
systems, the beacons can be broadcasted using the downlink 
channel. In TDD systems, the beacons could be given prece- 
dence over the user transmissions using the strategies similar 
to the ones based on SIFS/DIFS in IEEE 802.11 and related 
wireless standards, thus avoiding collisions. 

Finally, the proposed protocol could easily be extended for 
the scenarios with continuous arrivals, in the same way that 
related framed SA protocols operate jT], JU, iS). The protocol 
operation in this case would be divided into phases; all new 
arrivals occurring during a phase are backlogged, waiting for 
the next phase to start. In this case, provided that there is 
a correlation between the number arrivals over the phasefl 
the proposed estimating procedure could be straightforwardly 
modified to exploit it. We leave this for further work. 
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